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1. Introduction 

There has recently been a great deal of interest in investigating deformed thermodynamic 
systems at the classical level. Such deformed theories are believed to deal with the 
statistical behavior of complex systems, whose underlying dynamics is spanned in a 
multi-fractal phase-space, governed by long-range interaction and/or long-time memory 
effects [U [21 [3J H]. A possible fruitful mechanism capable of generating a deformed 
version of the classical statistical mechanics consists of replacing, in the Boltzmann- 
Gibbs distribution, the standard exponential with its deformed version, accordingly 
postulating a deformed entropic form which implies a generalized thermostatistics 
theory. In this manner, some noteworthy generalizations of the standard statistical 
mechanics have been proposed [21 EE IE] and their physical consequences are currently 
under investigation [HI [TO], [H] . 

On the other hand, quantum algebra and quantum groups have been the subject 
of intensive research in several physical fields from cosmic strings and black holes to 
the fractional quantum Hall effect and high-T c superconductors [12]. From the seminal 
work of Biedenharn [13] and Macfarlane |14j . it was clear that the g-calculus, originally 
introduced by Heine [15] and by Jackson [16] in the study of the basic hypergeometric 
series [T71 HE], plays a central role in the representation of the quantum groups with a 
deep physical meaning and not merely a mathematical exercise [T9l [20]. In this context, 
it was shown in Ref.[2T], that a natural realization of the thermostatistics of g-deformed 
bosons and fermions can be built on the formalism of g-calculus. 

Furthermore, it is remarkable to observe that the g-calculus is very well suited for 
to describe fractal and multifractal systems. As soon as the system exhibits a discrete- 
scale invariance, the natural tool is provided by Jackson g-derivative and g-integral, 
which constitute the natural generalization of the regular derivative and integral for 
discretely self-similar systems [22J. In fact, it was shown that g-integral is related to the 
free energy of spin systems on a hierarchical lattice [23] . 

In the recent past, some tentative to constructing a classical counterpart to the 
quantum group and g-deformed dynamics have been investigated [21]. Most recently, 
in [22], a g-deformed Poisson bracket has been developed whose underlying algebra 
arises from the quantum group theory appears to be invariant under the action of 
the g-symplectic group. This generalization implies a classical deformed dynamics and 
a deformed Fokker-Planck equation [26J whose stationary solution can be expressed 
in terms of the 6as«c-exponential, the g-analogue of the exponential function in the 
framework of the basic hypergeometric series (henceforward, we will use the term basic- 
exponential and 6asic-thermostatistics to avoid any confusion with the g-exponential 
customarily employed in the Tsallis' thermostatistics formulation [2]). 

The previous investigations raise the interesting question whether the g-calculus and 
the 6asic-exponential can be introduced, as a starting point, for the study of a deformed 
statistical thermodynamics at the classical level. Just as the quantum g-deformation 
plays a crucial role in the interpretation of several complex physical systems, we expect 
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that a classical g-deformation of the thermostatistics can be relevant in several physical 
applications. It is thus worthwhile to investigate the structure of a classical statistical 
mechanics where the probability distribution function is given by employing the basic- 
exponential. This investigation represents the primary goal of this paper. 

Relating to the previously quoted generalized statistical theories existing in the 
literature [21 E] , it must be stressed that the deformed exponential functions are 
very different from the one we are introducing in the present work. The difference 
arises mainly in the asymptotic behavior of the frasic-exponential which fails to show a 
power-law tail. Nevertheless, some mathematical peculiarities exhibited by the basic- 
exponential makes it relevant in the construction of a generalized statistical mechanics. 
First among them, the existence of a natural cut-off in the energy (velocity) spectrum 
which appears to be closely related to the presence of long-range interactions among 
constituents of the system. To the best of our knowledge, this was first suggested in 
Ref.[27j in order to overcome the infinity arising in the mass and radius of an isothermal 
globular cluster. More recently, concerning several physical systems, there have been 
many investigations [281 [29], EH! [311 [321 [331 Ell E2] on the relevance of the presence of 
an energy cut-off behavior in the particle distribution functions. 

Finally, we should stress that the fraszc-exponential introduced in the present work is 
also employed in the formulation of other interesting generalizations of physical theories 
such as, for instance, in the g-deformed Schrodinger theory describing the deformed 
quantum harmonic oscillator (36] or in the g-deformed theory of quantum coherence of 
bosons [33 EE]. 

Our paper is organized as follows. In Section II, for convenience of the reader, 
we review the main mathematical properties concerning the quantum algebra of real 
numbers. We introduce the frasic-exponential by means of a power series and derive its 
main proprieties which will be used in the formulation of the theory. In Section III, we 
introduce the basic-entropy and, by means of a g-version of the variational principle, 
we determine the deformed distribution, for each case in turn, of a microcanonical, 
canonical and grand canonical system. The thermodynamic structure of the theory is 
explored in Section IV, while in Section V, we apply the present formalism to a system 
of g-interacting particles whose q — > 1 limit reduces to the free system. Conclusions are 
reported in Section VI and we end the paper with two mathematical Appendices. 

2. Mathematical background 

We shall begin by recalling the main features of the g-calculus for real numbers. It is 
based on the following g-commutative relation among the operators x and d x , 



with q a real and positive parameter. 

A realization of the above algebra in terms of ordinary real numbers can be accomplished 



d x x = 1 + q x <9 : 



X 1 
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by the replacement 

x -> x , (2.2) 
d x ^V Xy (2.3) 

where V r is the Jackson derivative 1461 defined as 



„xd x _ 1 

= T r- • (2.4) 

X {q-l)x ^ 1 

Its action on an arbitrary real function f{x) is given by 

V.f(*) = fl f- flx) . (2-5) 
[q — 1) x 

The Jackson derivative satisfies some simple proprieties which will be useful in the 
following. For instance, its action on a monomial f(x) = x n is given by 

V x x n = [n}x n - 1 , (2.6) 

and 

V x x- n = - [ ^-t— , (2.7) 

q n x n+l V J 

where n > and 

q n - 1 

[»] = ^ , (2.8) 

are 6aszc-numbers. By linearity, we can extend the action of Jackson derivative to 
a generic polynomial. Moreover, we can easily verify the following g-version of the 
Leibnitz rule 

V x (f(x)g(x))=V x f(x)g(x) + f(qx)V x g(x) , 

= V x f(x)g(qx) + f(x)V x g(x) . (2.9) 

A relevant role in the g-algebra, as developed by Jackson, is given by the basic- 
binomial series defined by 

(x + y)^ = (x + y) (x + qy) (x + q 2 y) . . . (x + q n ~ l y) 

n 

= E [r] q r(r - 1)/2 x n - r y r , (2.10) 



r=0 



where 



- W (2 11) 

[r]![n-r]! ' 1 ' 

is known as the g-binomial coefficient which reduces to the ordinary binomial coefficient 

in the q 1 limit [18J. We should remark that Eq. (12.11j) holds for < r < n, while it 

is assumed to vanish otherwise and we have defined [n]\ = [n] [n — 1] . . . [1]. Equation 

fl2.10p can be easily generalized to an arbitrary polynomial as shown in Appendix B. 

Remarkably, a g-analogue of the Taylor expansion has been introduced in Ref. [16] by 

means of a basic-binomial (12.101) as 

(x-o)^ (x-o)W 

/w=/w+ L 1 ^i'./wL+ L 1 ^^/wL+... , (2.i2) 
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where V 2 = V x V x and so on. 

Consistently with the g-calculus, we also introduce the frctszc-integration 

A ° oo 

/ f(x)d q X=J2\Kf(K) , (2.13) 
n=0 

where A q X n = A n — A n+ i and A n = A q n for < q < 1 whilst A g A n = A n _i — A n and 
A n = \ q n l for g > 1 [T71 [18j [22]. Clearly, Eq. fl2.13l) is reminiscent of the Riemann 
quadrature formula performed now in a g-nonuniform hierarchical lattice with a variable 
step AqA n . It is trivial to verify that 



X 



-DxJ f(y)d q y = f(x) , (2.14) 
o 

for any q > 0. 

Let us now introduce the following g-deformed function defined by the series 

which will play the main role in the framework we are introducing. The function 
( I2.15P defines the fraszc-exponential, well known in the literature since a long time ago, 
originally introduced in the study of basic hypergeometric series [I5j [TFl [18]. In this 
context, let us observe that definition (I2.15P is fully consistent with its Taylor expansion, 
as given by Eq. ( 12.121) . 

In the following, we briefly review the main algebraic properties of the fraszc-exponential 
useful for our developments. 

The feasic-exponential is a monotonically increasing function, dE q (x)/dx > 0, convex, 
d 2 E q (x) / dx 2 > 0, with E 9 (0) = 1 and reducing to the ordinary exponential in the q — > 1 
limit: Ei(x) = exp(x). 

An important property satisfied by the g-exponential can be written formally as [18] 



E q (x + y) = E q (x)E 1/q {y) , (2.16) 

where the left hand side of Eq. ( 12. 16f) must be considered by means of its series expansion 
in terms of fraszc-binomials: 

fc=o 

By observing that (x — x)^ = for any k > 0, since (x — x)^ = 1, from Eq. (l2.16l) we 
obtain 

E q (x)E 1/q (-x) = 1 . (2.18) 

From the above relations, we easily deduce that, if q < 1, the 
series (12.151) converges for all finite values of x < 1/(1 — g), otherwise, if g > 1 the 
series converges for x > q/(l — q). 
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Thus, we can summarize the asymptotic behavior of the fraszc-exponential as 

1 

1-q 



E,(-oo) = , E, : = +oo , if q < 1 , (2.19) 



E q [y^j) = ' E <?(+°°) = +°° > if Q > 1 , (2.20) 

where it is important to mention that the first expression of Eq. (I2.20p defines a cut-off 
condition of the fraszc-exponential in the region relevant for our following developments. 
Among many properties, it is important to recall the following relation [18] 

V x E q (ax) = aE q (ax) , (2.21) 

and its dual 

E q (ay)d q y=-[E q (ax)-l]. (2.22) 
o 

It must be pointed out that Eqs. (I2.2ip and (12.221) are two important properties of the 
fraszc-exponential which turns out to be not true if we employ the ordinary derivative 
or integral. In particular, from Eq.f l2.21j ) for a = 1, we can obtain the further useful 
relation 

E q (qx) = [l + (q- l)x] E q (x) . (2.23) 
Moreover, from Eq. ( 12.221 ) with a < 0, we have 

E q (ay)d q y = -- , (2.24) 
a 
o 

where x max — > +oo in the q < 1 case, while x max = qj (q—1) in the q > 1 case, accounting 
for the cut-off condition (I2.20p . 

In addition to the fraszc-exponential, we can also introduce its inverse function, the 
ktsic-logarithm Ln q (x), such that 

E q (hn q (x)) = Ln g (E,(a;)) = x , (2.25) 

which certainly exists because E q (x) is a strictly monotonic function. Many proprieties of 
the fraszc-logarithm follow directly from the corresponding ones of the 6aszc-exponential. 
For instance, Ln q (x) is a monotonic, increasing and concave function (dLn q (x)/dx > 0, 
d 2 Ln ? (x) jd x 2 < 0), normalized in Ln 9 (l) = and the asymptotic behavior given by 

Ln q {0) = -oo , Ln q (+oo) = —— , if q < 1 , (2.26) 

1-q 

Ln„ (0) = , Ln„(+oo) = +oo , if q > 1 . (2.27) 

1 — q 

Although a definition of Ln q (x) through a series is possible, it appears to be a nontrivial 
task to write it in an easy form and, to the best of our knowledge, there are no definitive 
results in the literature (see for instance Ref.jH]). 

We conclude this Section by remarking that alternative definitions of the basic- 
exponential, by means of a different definition of 6asic-numbers, has been widely 
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employed in the literature [42 J . Among the many, we may point out the choice based on 
the symmetric definition [n] s = {q n —q~ n ) / (q—q^ 1 ). The corresponding symmetric basic- 
exponential may be defined on the whole real region (—00, +00) and has the symmetry 
q — > 1/q. As a consequence, the symmetric frasze-exponential does not present the cut-off 
feature. 



3. Basic-entropy and its distribution 

Equipped with the fraszc-functions, our aim is to formulate a statistical mechanics based 
on the formalism of the g-algebra and to study its main physical implications. 
On the basis of the above mathematical framework, it appears natural to generalize the 
Boltzmann entropy to the following form 

S q (p) = ~ Jp(X)Ln q p(X)d q X, (3.1) 

M 

where p(X) is probability distribution function labeled by a set of parameters A running 
on the manifold Ai, eventually identified with the phase space, which define the 
accessible states of the system. Henceforward we adopt units where the Boltzmann 
constant k-Q = 1. 

Equation (13.11) resembles the well-known Boltzmann- Gibbs entropy S BG (p) through the 
replacement of the logarithm with the basic- logarithm. Clearly, the function S q (p) 
reduces to the standard entropy S BG (p) = Si(p) in the q — > 1 limit. In the following, 
we shall refer to the function ( 13. ip as 6aszc-entropy. 

A way to obtain the equilibrium distribution from the entropy S q (p), given a set of 
M + l constraints &j(p) with j = 0, . . . , M, can be accomplished through the following 
variational problem 

5T{p) = , (3.2) 

with 

^(P)= fs fl (p)-EA*i*i(p)) > ( 3 - 3 ) 

where /1* are the Lagrangian multipliers associated with the constraints 
Quite generally, such constraints can be written in 



<t> j (p) = J ( j> J (X)p(X)d q X, (3.4) 

M 

representing the mean value of the quantities 4>j(X) which are identifiable with suitably 
physical observable. In particular, for </>o(A) = 1, Eq. (13.40 give the normalization of the 
distribution function. 

In order to handle easily the variation problem of Eq.( 13.2p . let us introduce the ansatz 

p(X) = E q (-f(X)) . (3.5) 

We recall that E g (— x) is a strictly monotonically decreasing function. This means that 
it reaches its minimum at points which maximize the functions /(A). 
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Actually, the problem Eq. fj3.2j) can be replaced by the following equivalent q- variational 
problem 

6 9 f(f)=V f f(f)A q f = 0, (3.6) 
with A g f — (q — l)f and 

Hf) = Hp(f)) = I E g ( - /(A)) (/(A) - J> M*)) ^ • (3-7) 



M 3 

It is show in Appendix A that, according to the g-algebra, both the calculus (I3.2p and 
(13. 6j) give substantially the same result (apart a redefinition of the Lagrange multipliers, 
fij — > fi*, which, accounting for the relevant constraint equations, has no any effect on 
the expression of the final distribution). Equation (13. 6ft can solve more speedily as 



V f jE q (- /(A)) (/(A) - 5>i (pj(\fj d q \ 
. M i 

/E 9 ( - /(A)) (qf(X) - 1 -£/*,■ &(A)) d 9 A 



A J = , (3.8) 



which implies the following relation 

/(A) = g- 1 (l + 5:/i,^ 

and accounting for Eq. (13.51) we obtain the general solution in the form 

p(X) = E q (^-q- 1 (l + E^^iW)) • 



(3.9) 



(3.10) 



It is easy to verify that this expression reduces, in the g — ^ 1 limit, to the ordinary Gibbs 
distribution. Let us consider separately three main cases: the microcanonical system, 
the canonical system and the grand canonical system. 

3.1. Microcanonical system 

We consider a closed system with a given fixed energy U, volume V and particle number 
N. In this case the system is only forced by the constraint 



p(X) d q \ 



(3.1.1) 



M 



which assures the normalization of the distribution. 

Before to proceed on, let us spend a few word about Eq. (l3.1.1l) to better understand the 
underlying physical meaning of the g-calculus. By taking into account the definition of 
the 6aszc-integral (I2.13p . for < q < 1, we obtain 



p(A) d q X = A (1 - q) p(A ) + qp(Xi) + q 2 p{X 2 ) + 



(3.1.2) 



M 



where Ao is a constant and Aj = q % X§. This expression shows two important features of 
the theory we are developing. First, the g-integral plunges, in a natural way, to consider a 
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g-deformed lattice whose amplitude of the elementary cell A ? A = A (1 — q) is shrunken, 
step by step, by the quantity q n . This is a substantial different situation with respect 
to the standard case obtained in the q — > 1 limit where any probability (differential 
probability) is multiplied by an equal (infinitesimal) quantities: dp(x) = p(x) dx. 
Second, the parameter space assumes a fractal structure given by the rule X n = q n Ao- 
We recognize a self-similarity in the parameter space since, starting from any level N > 1 
the same structure X^ +n = q n Xn is discovered. 

Remarkably, by introducing a set of discrete probability distributions p n , related to p(X), 
as 

p n = A q X n p(X n ) , (3.1.3) 
equation (13.1. ID becomes 

oo 

£ft = l. (3-1.4) 

i=0 

Clearly, the same considerations hold also in the q > 1 case. 

In order to derive the expression forp(A) we introduce the following constrained entropic 
functional 

F(P) = ~J PW [Ln 9 p(A) + a] d q X , (3.1.5) 

M 

where a is the Lagrange multiplier associated with Eq. (13. 1.1 1) . 
By inserting the ansatz (I3.5P in Eq. (13. 1.51) we obtain 



Hf) = J E ? ( - /(A) J (/(A) - a) d q X . (3.1.6) 

M 

By evaluating the equation 5 q F(f) = 0, accounting for Eqs. (12.91) and (I2.2ip . we obtain 

5 q F{f) = E q ( - /(A)) (l + a -qf(X)) = . (3.1.7) 

The above expression implies the following equation 

Ln q p(X) + q- 1 (l+a) = , (3.1.8) 

which gives the microcanonical distribution in the form 

p(X) = E q (-q- 1 (l + a)) . (3.1.9) 

[We refer to Appendix B for a systematic derivation of Eq. (I3.1.9P . 

Since the above expression does not depend on A, by accounting for the condition ( 13. 1.1ft . 

we obtain the expected microcanonical uniform distribution 

P=^, (3.1.10) 
where the number of accessible states 

W(U, V, N) = / d q X , (3.1.11) 



M 



is related to the Lagrange multiplier through the relation a = —1 — qLn q (l/W), which 
is a function of the energy, the volume and the total number of particles in the system. 
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By following the steps described in the previous Subsection, we can derive the canonical 
distribution for an open system that exchanges energy with the surrounding. In this 
case, let us pose e = e(A) the energy of the system still labeled through A so that, among 
to the constraint (13.1.11) . we impose the further condition on the mean energy 



<A)p(e(A)KA = <6) 



(3.2.1) 



M 



(3.2.2) 



Accounting for the definition (12.131) this last condition becomes (0 < q < 1) 

A (1 - q) [e p(e ) + <?eip(e 1 ) + q 2 e 2 p(e 2 ) H ] = (e) , 

where ej = e(Aj) so that the g-calculus implies a fractal structure in the energy spectrum. 
With the position ( 13.1 .3D we realize that the constraint ( 13.2.11) on the mean energy is 
equivalent to the standard definition 



8=0 



(e) 



(3.2.3) 



although the fractal structure in the energy spectrum still holds being implicitly 
contained in the definition of the discrete probabilities Pi. 
After introducing the constrained entropic functional 



Hp) = ~ p(*W) Ln q p(e(X)) + a + (3e(X) 



d q X 



(3.2.4) 



M 



with a and (3 the Lagrange multiplier associated with Eqs.( l3.1.TT) and (13. 2.1ft . by 
imposing the ansatz (13.51) we are dealing with the variational problem 5 q J r (f) = 
where 



Hf) = /E,( - /(6(A))) [/(e(A)) - a -0e{\) 



dg\ 



(3.2.5) 



M 



and whose solution reads 

Ln q p(e) + q~ 1 (1 + a + (3 e 



. (3.2.6) 

According to the g-algebra described in Eq. (1B. 41) . with x = Ln q p(e), y = q^ 1 (1 + a) 
and z = q~ l (3 e, from Eq. (13.2.61) we derive the following canonical distribution 



P(e) = E 9 



-q'^l + a] 



E q (-/3 q e), (3.2.7) 

where (3 q = q~ l (3 (see Appendix B). By imposing the normalization condition ( 13.1.11) 
on the distribution p(e), we obtain 



p (e) = ^E q (-(3 q 



where Z q is the canonical partition function defined by 



Z q = E 1/q (q- 1 (1 +a))=J E,( - (3 q e(X)) d q X 

M 



(3.2.8) 



(3.2.9) 
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Trivially, Eg. (I3.2.8|) reduces to the canonical Gibbs distribution in the q — > 1 limit. 
We remark that, for q < 1 all the energy levels e can be occupied. For q > 1, however, the 
distribution (13.2.81) shows a cut-off in the energy spectrum due to the finite convergence 
radius of the function E q (x). This is an important consequence of the theory under 
investigation which limits the number of states accessible to the system. Its origin can 
be related, as generally accepted, to the presence of interactions among the parts of 
the system, whose effect is to reduce the volume of the phase space. This peculiarity 
is encountered also in other statistics mechanical models based on generalized entropic 
forms. 

In particular, Eq. (12.201) imposes the following limiting condition on the energy levels 



e < e m ax, where 



'' (3.2.10) 



mdX (? - 1) q 

Physically, this means that all the microscopic configurations of the phase space 
corresponding to an energy e beyond e max are statistically unattainable. We remark, that 
e max is a function of f3 q which plays the role of the inverse of a pseudo temperature. Its 
value is determined through Eq. (13.2.11) and it is expected, like in the undeformed theory, 
that (3 q decreases as the mean energy increases. It means that the cut-off condition plays 
a relevant role in those small systems whose mean energy is small compared with the 
typical energy values of the macroscopic systems. This is reminiscent of a quantum 
scenario although our system is a classical one. In Figure 1, we illustrate the behavior 
of the basic Boltzmann factor E q (—j3 q e) for different values of the deformation parameter 
q compared with the classical one given by q = 1. It is observed that for q < 1, high 
energy events are enhanced with respect to the standard case while, for q > 1, events 
are more and more inhibited when energy increases until it reaches the cut-off point, 
where p(e) = 0. 

In Figure 2, we show the frasic-entropy for a system of two levels with probabilities p 
and 1 — p, respectively, for different values of the deformation parameter. The dashed 
line represents the asymptotic curve reached for q ^> 1. 

Again, we stress that such a feature is also shown by other distributions obtained from 
physically motivated generalizations of the Boltzmann- Gibbs entropy. For instance, it 
can be shown [TO] that for suitably chosen deformation parameters, generalized entropies 
belonging to the two-parameter family of Sharma-Mittal [43J , which include among the 
others also the well known Tsallis' entropy [5] and the Renyi's entropy [44], generate 
probability distribution functions which exhibit a cut-off in their tails. Nevertheless, it 
is worthwhile to observe that almost all the members belonging to the Sharma-Mittal 
family, with some exception like the Boltzmann-Gibbs entropy and the Gaussian entropy, 
have an asymptotic power law behavior. This differs substantially from the asymptotic 
behavior shown by 6aszc-distribution, which is more similar to that of the stretched 
exponential. In this respect, the theory under investigation is not an alternative but 
complementary to the already existing generalized version of the statistical mechanics, 
since it can be relevant in the study of those complex systems which are not characterized 
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Figure 1. Generalized Boltzmann distribution E g (— /3 q e) for different values of the 
deformed parameter q. The solid line coincides with the standard exponential function 
exp(-/3e). 



by an asymptotic free-scale behavior. 



3.3. Grand canonical system 



Finally, let us investigate the grand canonical distribution describing an open system 
where energy and particles can be exchanged with the surrounding. This can be 
accomplished by imposing the following constraints 



Wp„(e(A)KA = l, (3.3.1) 
N M 

E fe(X)p If (e(X))d q X=(e) , (3.3.2) 

N M 

J2[Np N (e(\))d q \=(N), (3.3.3) 

AT J 



N M 



on the the normalization, the mean energy and the mean particle number, where 
N — 1, . . . , oo enumerate the particles contained in the system. 

The distribution function is found just as in the former cases. We form the constrained 
entropy 

HP) = "E / P JV (e(A))[Ln (? p iV (e(A)) + a + (3e{\) + 1 n] d q \ , (3.3.4) 

N M 
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Figure 2. Plot of basic-entropy, for different values of the deformation parameter, for 
a system with two levels. The dashed line represents the asymptotic curve reached for 
g> 1. 



where a, (3 and 7 are the Lagrange multipliers associated with the constraints (13.3.11) - 
(13.3.31) . respectively. By evaluating the equation 8 q T{f) = 0, with 

Hf) = E [M ~ [fMV) - a -/3e(A) - 7 JV] <**A , (3.3.5) 

N M 

obtained from Eq. (l3.3.4p by using the ansatz p N (e) = E 9 ( — f N (efj, we derive the 
following result 

Ln g p N (e) + q~ 1 (l + a + Pe- fxPN) = , (3.3.6) 

where we set 7 = — According to the g-algebra (IB. 41) with x = Ln q p N (e), 
y = q~ l (1 + a), z = q~ l (3 e and u = —q^ 1 n(3 N, from Eq. ( 13.3.61) we obtain the grand 
canonical distribution function in the form 

p N (e) = ^E q (-(3 q e)E q ( f if3 q N) , (3.3.7) 

where 

Z q = E 1/q (q-\l + «))=£/ E ff ( - P g e(X)) E^faN) d q \ , (3.3.8) 



N M 



is the grand partition function. 
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4. Some thermodynamic relations 

In the following we shall investigate the thermodynamic structure of the theory. It 
is shown that some basic relations of the standard theory can be transcribed in a 
straightforward manner in the present formalism establishing, in this way, the ground 
for a generalized classical thermodynamics based on the framework of the g-algebra. 
To begin with, we observe that by multiplying Eq. (I3.3.6|) by p N (e) and taking into 
account all the constraints imposed on the system, we obtain the relation 

S q ({e), (N),V) = Ln 1/q Z q + /3 q (e)-til3 q (N) , (4.1) 

where we set Lni/ 3 Z q = q^ 1 (1 + a) as it follows from Eq. (I3.3.8j) . Equation (14. ip mimics 
the standard relation S — In Z + j3U — fi f3 N which is recovered in the q — > 1 limit. 
If one is willing to identify the quantity (3 q with the inverse of the temperature (3 q = 1/T, 
Eq. fl4.ll) can be rewritten in 

(e)=TS q -TLn 1/q Z q + fx (N) . (4.2) 

On the other hand, in analogy with standard thermodynamics, we can introduce the 
pressure P and the volume V of the system by requiring that all the thermodynamics 
variables are functionally related through the following relationship 

(e)=TS q -PV + fi(N) , (4.3) 

which is a constitutive equation for the theory under investigation and can be identified 
with the g-analogue of the Euler's equation [35]. By comparing Eqs. fl4.2p and (14. 3p we 
are thus encouraged to define the pressure through the relation 

PV = TLn 1/q Z q , (4.4) 

which represents the g-generalized state equation for a system described by the basic- 
entropy. 

The partition function f)3.2.9p . as well as the grand partition function (13.3.81) . are 
useful tools to evaluate the statistical proprieties of the system. In fact, by evaluating 
the Jackson derivative of Z q as 

V Pq Z q = V Pq J E, ( - (3 q e(A)) d q \ = - J e(A) E, ( - (3 q e(A)) d q X , (4.5) 

M M 

where Eq. (12.211) has been used, we obtain the result 

{e) = ~V Pt Z q , (4.6) 

which, in the q — > 1 limit reduces to the well-known relation (e) = —{dZ/df3)/Z. 
Similar results can be obtained starting from the function Z q as follows 

(e) = -^V Pq Z q , (4.7) 
(N) = -±-V IM Z q . (4.8) 

PqZ q 
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Incidently, by using the prescription 

e^e + 5A(e), (4.9) 

where A(e) is an arbitrary physical observable, the expectation value of A can be 
obtained in 



(4.10) 

5=0 

In this sense the partition function encode all the statistical information contained in 
the system. 

Another physically relevant quantity is obtained starting from the expression 

^<c> = ^(^r/e(A)E B (-/9,c(A))^A] , (4.11) 



q 'M 



where we have explicitly indicated the dependence on f3 q in the partition function. By 
taking into account Eq. fl2.23p . we can derive the result 



V ^ = - !-(?-!) t W ^ • (4 ' 12) 
where (e 2 ) = / e 2 p(e) d q X. 

In the q — > 1 limit, we obtain the classical relation d (e)/d (3 = — o~ 2 . 
The quantity a q ^ £ measures the fluctuation of the energy of the system around its mean 
value (e). It is observed that, compared to the classical case, such fluctuations are 
reduced for q < 1 and are enlarged for q > 1. 

Alternatively, by mimicking the classical definition, one can introduce the heat capacity 
by 

C v = -(3 2 q Vp q (e) = (P q a q , £ ) 2 , (4.13) 
so that, the relative width of the fluctuations in energy is given by 

7T = TT~\^ ■ ( 4 - 14 ) 

Similar derivations can also be made for the grand canonical case where, together with 
Eq. fl4.12l) we can obtain also the further relation 



V»(N) = P q ,JP u {N « 2 /m = <n , (4-15) 
l + (q-l)fi/3 q (N) 

with (iV 2 ) = J2n I N 2 p N (e) d q X, which measures the deviation from the mean value (N). 



5. Basic-ideal gas 

In this Section, to illustrate the consequences of our generalized thermodynamical model, 
we discuss a simple 6as«c-noninteracting particle gas which reduces, in the q — > 1 limit, 
to the well known ideal gas. 
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Let us start from the following Hamiltonian H describing a system of N identical 
particles 

N 2 

tfo(p) = E^T> (5-1) 

where p = (p l3 . . . , p N ) is the 3 N- vector momenta. In order to be consistent with 
our model, we require that the momenta p t obey the g-algebra discussed in Appendix 
B. In this case, we can verify that the canonical distribution f(p) obtained from the 
6aszc-entropy (13.11) with the mean energy constraint 

j H{p)f{p)dfxdl N P =(e) , (5.2) 

M 

can be factorized as follows 

N 
i=l 

Here, /(pj is the probability distribution function of a single particle 

/<K> = ^E.(-Al?L) , (5.4) 
with p 2 = p 2 ■ + p 2 ■ + p 2 ■ and 



2 \ /o \ 3/2 



Z q>i = I E q \ —p q 1 djd djx = V lA \ ~3 ) » (5-5) 



M 



in units where Ti = 1. In the above equation, we have posed V^ q = j M d? q Xi the fractal 
volume occupied by the z-th particle and m q = mA q where 

1 °° 

A q = —J E q (-x 2 ) d q x , (5.6) 



is a constant reducing to unity in the q — > 1 limit. The canonical partition function of 
the whole system is given by 

N f2nm V N/2 

^=n^ = ^(^J • (5-7) 

We may note that Eq. f!5.4p can be interpreted as the g-deformed version of the Maxwell- 
Boltzmann distribution. It differs formally from the well-known classical distribution by 
the mere replacement of the standard exponential with its g-deformed generalization, 
consistently with the g-algebra underlying the mathematical structure of the theory. 

Employing the distribution f(p), we can compute the mean value of any observable 
associated with the system. In particular the mean energy is given by 



L2 



From Eq. fl5.8p . we recognize in the q — > 1 limit the well-known 
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result (e) = |iVT. We observe that 

(e) -> +00 for q < , and (e) ->• for g > 1 , (5.9) 

which is a consequence of the dependence of f(p) on g. In fact, for smaller and smaller 
q the tail of the distribution is enhanced so that particles with high energy give a larger 
contribution. In contrast, for larger and larger q the cut-off inhibits the occupation by 
the system of the space-phase cells with high energy and only particles with lower and 
lower energy contribute to (e). In particular, for q ^> 1 the distribution f(p) — > 5(0) 
so that only the fundamental level is populated. Similar arguments can be applied to 
justify the expression of the heat capacity given by 



(5.10) 



which is a constant as in the undeformed classical case, but it is a monotonically 
decreasing function reducing to zero for q ^> 1. 

It is important to outline that the free particle gas with the Hamiltonian (15.11) has 
the g-deformed particle distribution (15.41) only if we require, as a crucial assumption, 
that the particle momenta obey the g-algebra. On the other hand, if we do not employ 
the appropriate g-algebra, the same free g-deformed particle distribution (15.41) can be 
obtained by assuming the following iV-body interacting Hamiltonian 



This suggests that the effects of the basic- deformation of a /ree-ideal gas can be viewed 
as an effective interaction described by the Hamiltonian ( 15. lip . We remark that in the 
q — > 1 limit, the g-algebra reduces to the ordinary algebra used in the Hamiltonian 
(15. lip . Furthermore, the same Hamiltonian (15. lip reduces to the Hamiltonian of a 
iV-free particles system. 

It may be important to clarify this point. In several papers [8], 061 HH 08], it has 
been suggested that, starting from a deformed exponential derived through physically 
and/or mathematically justified arguments, it is possible to introduce a deformed sum 
in order to reproduce, in a deformed fashion, the well-known multiplicative rule of the 
standard exponential exp(x + y) — exp(x) exp(y). It has been conjectured (see for 
instance [50]) that such a deformed sum can be employed, on physical grounds, 
to take into account the complex interactions arising among the many-body colliding 
particles of a nonlinear medium. This has been discussed, for instance, explicitly in the 
Tsallis-entropy framework [51J . In that case in fact, for the deformed sum of the energy 
values £ A and £ B belonging to two different subsystems A and B, it has been assumed 
that the expression describes the g-sum 

£ A ® g S B = S A + S B + ^ S A S B , (5.12) 

P 

and correspondingly the g-Boltzmann factor factorizes according to 

exp q (-[3(S A ® q S B )) = ex P(? ( - (3S A ) ex P(? ( - (3 S B ) . (5.13) 
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The same situation occurs in the Kaniadakis-entropy framework [7] where, by assuming 

£ A © £ B = £ A sjl + (k(3£*) 2 + £ B \Jl + (k(3£ a ) 2 , (5.14) 

for the energy levels, it has been shown that the K-Boltzmann factor factorizes in 

exp {K} ( - (£ A © £ B )) = exp {K} ( - PS A ) exp {K} ( - $ £ B ) . (5.15) 

Of course, all of this can be also reproduced within the formalism employed in the 
present work. In fact, we can verify that the following deformed sum 

£ A ®£ B = -^Ln q (E 1/q (p£ A )E 1/q (p£ A )) , (5.16) 

which reduces to the ordinary sum in the q — > 1 limit, fulfils the factorization rule 

E q (-f3(£ A ®£ B )) = E q (-f3£ A )E g (-(3£ B ) , (5.17) 

for the ftaszc-Boltzmann factor. We easily recognize in the definition fl 5 . 1 6 [) the origin of 
the Hamiltonian (15.111) . Although it is beyond the scope of the present paper, it is also 
possible to show that the basic-sum (15. 16ft obeys all the axiomatic properties where a 
well defined sum must satisfy associativity, commutativity, existence and uniqueness of 
the inverse element and of the identity element. However, by applying the appropriate 
g-algebra introduced starting from the g-binomial expansion , the deformed sum x © y 
is replaced, in a natural way, by the ordinary sum x + y. This aspect makes the basic 
thermostatistics formalism very interesting, because of the structure of the considered 
deformation implies a close and consistent realization in the well known mathematical 
framework of g-calculus. 

6. Conclusion 

In this paper we have studied a possible generalization of the thermostatistics theory 
of a classical system based on the g-deformed algebra. Starting from the definition of 
the ftetszc-exponential, we have introduced a generalized entropic function and derived, 
by means of a consistent g-variational principle, a deformed probability distribution 
function which differs from the standard Gibbs distribution by the replacement of the 
ordinary exponential function with its generalization furnished by the 6aszc-exponential. 
We have performed a preliminary investigation of some fundamental thermodynamic 
relations which are preserved consistently with the formalism of the g-calculus. 

On physical grounds, it has been demonstrated that the distribution arising in this 
model exhibits a cut-off in the energy spectrum which is generally expected in those 
systems whose underlying dynamics is governed by long-range interactions. Such a 
feature has been observed also in other distributions proposed in the literature, obtained 
from generalized versions of the Boltzamnn-Gibbs entropy. What is different here, is 
the asymptotic behavior of the distribution derived in this paper which does not match 
with the power-law behavior typically shown by the other generalized distributions. 

We have studied, within the present formalism, a iV-body system of interacting 
particles described by the Hamiltonian (15.111) . whose interaction vanishes in the q — > 1 



Basic- deformed thermo statistics 



19 



limit. By construction, the canonical distribution function of this system is formally 
equivalent to the one derived starting from the Hamiltonian (I5.ip . describing a system of 
non-interacting particles, where the momenta obeying the g-algebra originates from the 
g-binomial expansion (12.101) . In this sense, the toy model described by the Hamiltonian 
(15. lip can be considered as the g-analog of the ideal gas. 

Being Jackson g-derivative and g-integral the natural tools for describing discrete- 
scale invariance [221 [22] , we expect that this study can be a very useful starting point 
on the feasic-thermostatistic theory which can be strictly related to critical phenomena 
(such as growth processes, rupture, earthquake, financial crashes) with the existence 
of log-periodic oscillations deriving from a partial breakdown of the continuous scale- 
invariance symmetry into a discrete-scale-invariance symmetry, as occurs for instance in 
hierarchical lattice [52j. In fact, as mentioned in the Subsections 3.1 and 3.2, our study 
on the feasic-thermostatistics incorporates implicitly a self-similarity in the parameter 
space, labeling the phase space of the system, and, consequently, a fractal structure 
in the energy spectrum emerges. In this context, it is remarkable to observe that, for 
example, the specific heat corresponding to systems with deterministic fractal energy 
is known to present log-periodic oscillations as a function of the temperature around a 
mean value given by a characteristic dimension of the energy spectrum [53l IMf [55] . A 
detailed study of these critical phenomena in our formalism lies out the scope of this 
paper and will be the matter of future investigations. 



APPENDIX A. 



As known, a possible way to obtain the stationary distribution of a system governed by 
a given entropy, under a set of suitable physically constraints, follows by means of the 
variational calculus on the constrained entropic form. In the case under investigation, 
accounting the g-algebra underling our formalism, we have to deal with the following 
problem 

5F(p) = 0, (A. 1) 

where, according to Eq. (l3.3p 



Hp) 



p(A) [Wp(A)) + X>ifc(A) 



d q X 



(A. 2) 



M 



Without loss of generality, we pose /ig = g~ 1 (l + /io) — 1 and /i* = q~ x fa for i — 1, ... M, 
where fij are the Lagrange multipliers of to the M + 1 constraints (13 Ah introduced in 
the Section 3. 

By using the ansatz (13. 5p in ^(p) = J-(p(f)), Eq.l lA. ll) can be computed as follows 

E, 



f(X)-th(X))d q X 



/(A)-5>*<MA)] E q (-f(X))d q X 



M 
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d 

~di 



M 



J [/(A) + 1 h(X) - £ /i* 0j(A)J E g ( - /(A) - 1 h(X)) d q X 

£E 1/q (-th(X) 



t=o 



/»(A) + (/(A)+f/»(A)-5>^(A) 



E v ,(-*/i(A) 



E ff (-/(A))d,A, 



t=o 



(A. 3) 



where E g (-/ -th) = E q (-f) E 1/q (-t h) according to Eq. fl2~T6l) . 
By taking into account the definition (I2.15P we have 



It El/ « { 



-MA) 



(A. 4) 



^(A)) = Eo T (-MA))r- 1 

n=l l' 4 J- t=0 

so that from Eq. ll A. 31) we obtain 

= / h(X)(l - f(X) + X>*<MA)) E,( - /(A)) d q X = , (A. 5) 



Accounting for the arbitrariness of the function h(X) this last equation implies 

/(A) = l + 5»i(A) , (A. 6) 

i 

in accordance with Eq. (13. 9p given in the text. 
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We present a generalization of the algebra (12.101) to a trinomial with the purpose of 
extending the factorization formula of the fraszc-exponential. 
First, let us briefly review the derivation of Eq. ( 12.161) . 

This can be shown easily by considering the Cauchy product among the series (12.151) 
and its analogue for q — ► 1/q. We obtain 



E q (x) E 1/q (y) 



1+1 — + -*- 
Ul]! [1]! 



X~ 



+ 



xy 



+ 



[2]! [1]![1]! [2]! 



+ . . . + 



+ 



x n x n-l y q X n - 2 y 2 ^ g n(n-l)/2 y n 



+ 



\[n}\ [n-l]![l]! [n-2]![2]i [n\\ j 

which, by means of Eq. fj2.10j) . can be rewritten in the form 

e, (*) E 1/t(y) = i + L^L. + L_|L_ + . . . + L_|_ + 

and coincides with the definition of E q (x + y) given in Eq. (12.17| ). 
To generalize this result, we consider the following g-binomial expansions 
m = 1 



[X + z) 
(x + z t 
(x + z) 
(x + z 



,(1) 

(2) 
,(3) 



X + Z 
.,2 



x 2 + [2] x z + q z 2 , 
x 3 + [3] x 2 z + q [3] x z 2 + q 3 z 3 



,(B. 1) 



(B. 2) 



(B. 3) 



Basic- deformed thermo statistics 21 

and so on. By redefining x — > x + y and consequently x n — > (x + y)^, we obtain 

((s + y)+z) (0) = l, 

{{x + y) + z) (1) = x + y + z , 

((x + y) + z) {2) = (x + y) {2) + [2] (x + y) {1) z + qz 2 

= x 2 + [2]xy + qy 2 + [2}x z + [2}y z + qz 2 , 

((ar + y) + z) (3) = (x + y) (3) + [3] (x + y) {2) z + q [3] (x + y) {1) z 2 + q 3 z 3 
= x 3 + [3] x 2 y + q [3] xy 2 + q 3 y 3 + [3] x 2 z 

+ [2] xy z + qy 2 z + q[3] x z 2 + q[3}yz 2 + q 3 z 3 , (B. 4) 

which implies the following factorization rule for the feaszc-exponential 

E g (x + y + z) = E q (x + y) E 1/q (z) = E q (x) E 1/q (y) E 1/q {z) . (B. 5) 

On the other hand, starting from the g-binomial expansion ( IB. 31) and by posing 
z n — > (y + zfy q 0we form the g-trinomial expansion (x + (y + z)i/ q ) <yn " > which implies the 
following decomposition 

E q {x + y + z) = E q {x) E l/q (y + z) = E q (x) E l/q (y) E q (z) . (B. 6) 

Other possible factorization rules can be realized through the introduction of suitable 
g-deformed algebras, as can be seen by inspection. Extension to an arbitrary polynomial 
can be also easily obtained. 

In the following, let us show the use of the above algebra in the derivation of the 
distributions ( 13. 1.9ft . ( 13. 2 .TP and (13.3.71) . Starting from the equality 

x + y = , (B. 7) 

and employing the g-algebra (IB. 3j) we can construct the g-binomial (x + y)^ = which 
of course vanishes for all n > 0. Consequently, according to the definition ( 12. 15ft . by 
dividing Eq. (1B. 7j) by [n]\ and summing up over n = 0, . . . , oo, we obtain 

E q {x + y) = E q {x)E 1/q {y) = l . (B. 8) 

In particular, by applying this result to Eq. (l3.1.8l) which has the form ( IB. 71) . with 
x = Ln g p(e) and y = g _1 (l + a), we obtain 

E g (Ln g p(e)) E^q' 1 (1 + a)) = 1 , (B. 9) 

so that 

p(e)=E q (-q- 1 (l + a)) , (B. 10) 

and according to the property ( 12.181) we obtain Eq.( l3.1.91) given in Section III- A. 
In the same manner, from the equality 

Ln q p(e) + q- 1 (l + a + (3e)=0 , (B. 11) 

| We have introduced the index 1/q to indicate the replacement q — > 1/q in the expansion of the 
g-binomial given in Eq. (I2.10p 
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given in Eq. (I3.2.6j) . employing the g-algebra (IB. 4p with x = Ln q p(e), y = q 1 (1 + a) 
and z = q~ l (3 e, it follows that 

E ? (Ln g p(e))E 1/? ( g - 1 (l + Q; ))E 1/? (g- 1 /5e) = 1 , (B. 12) 

or equivalently 

ft = E 9 (-g- 1 (l + a))E g (-g- 1 /3e) , (B. 13) 

and by using again the property (12. 18ft it can be written in the form ( 13.2.71) given in 
Section III-B. Similar arguments can be applies to obtain the distribution (13.3.71) by 
employing the appropriate g-algebra with x = Ln q p(e), y = q^ 1 (1 + a), z = q~ x (3 e and 
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